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Abstract: The amoeba of a complex hypersurface is its image under a logarith- 
mic projection. A number of properties of algebraic hypersurface amoebas are 
carried over to the case of transcendental hypersurfaces. We demonstrate the po- 
tential that amoebas can bring into statistical physics by considering the problem 
of energy distribution in a quantum thermodynamic ensemble. The spectrum 
{sk} C Z" of the ensemble is assumed to be multidimensional; this leads us 
to the notions of a multidimensional temperature and a vector of differential 
thermodynamic forms. Strictly speaking, in the paper we develop the multidi- 
mensional Darwin and Fowler method and give the description of the domain of 
admissible average values of energy for which the thermodynamic limit exists. 



1. Introduction 

The amoeba of a complex hypersurface V defined in a Reinhardt domain is the 
image of V in the logarithmic scale. The notion of the amoeba of an algebraic 
hypersurface, introduced in [I], plays the fundamental role in the study of zero 
distributions of polynomials in C". Over the last decade, the amoebas proved 
to be a useful tool and a convenient language in the diverse questions; such as 
the classification of topological types of Harnack curves [2] , description of phase 
diagrams of dimer models J3j[4] , study of the asymptotic behavior of solutions to 
multidimensional difference equations [5]. Adelic (non-archimedeam) amoebas 
turned out to be helpful in the computation of nonexpansive sets for dynamical 
systems |B]. 

The main purpose of the present paper is to demonstrate the advantages of 
using the amoebas in statistical physics. As an example of such usage we consider 
the statistical problem of finding the preferred states of the thermodynamic 
ensemble when its spectrum is discrete. 
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In the classical formulation of this problem, which was studied by Maxwell, 
Boltzmann and Gibbs, it is assumed that the energy levels occupied by the en- 
semble systems form a discrete one- dimensional spectrum {£&} CN = {0, 1,2,.. .} 
(see, for example, 0H])- By contrast, we consider the case of the multidimen- 
sional spectrum {e^} C Z n , n > 1, and then such important notions as the 
temperature and the differential thermodynamic form become vector quantities. 

In fact, the major part of the paper is devoted to the generalization of the 
asymptotic Darwin- Fowler method [9,110], that gives a way to describe the state 
of a quantum thermodynamic ensemble with the multidimensional spectrum. 
For this purpose, we introduce the notion of an amoeba of a general (not only 
algebraic) complex hypersurface and describe the structure of the amoeba com- 
plement (Theorem [I]) . Next, we prove an asymptotic formula (Theorem [2| for 
the diagonal Laurent coefficient of a meromorphic function; the polar hypersur- 
face, its amoeba and the logarithmic Gauss mapping are significantly used in 
the proof. 

There are two main reasons motivating to apply methods of the theory of 
amoebas to the asymptotic investigation of the Laurent coefficient of a mero- 
morphic function in several complex variables. First, the connected components 
of the amoeba complement are in one-to-one correspondence with the Laurent 
expansions of a meromorphic function centered at the origin; and, moreover, de- 
fine their domain of convergence. Second, by the multidimensional residues, the 
asymptotics of the Laurent coefficient is represented by the oscillating integral 
over a chain on a polar hypersurface V. In the logarithmic scale, the critical 
points of the phase function of such integral comprise the contour of the amoeba 
of V. 

Thus, our generalization of the Darwin-Fowler method (Sect. [7]) is grounded 
on Theorems [l] and [2j Theorem [3] provides the asymptotics of the average val- 
ues for occupation numbers of energy Sk from a given spectrum. These average 
values are expressed by the Laurent coefficients of the meromorphic function 
constructed by means of the partition function of an ensemble. Although The- 
orem [3] requires tricky integration techniques, its statement is a quite expected 
generalization of the Darwin-Fowler results. This is not the case with Theo- 
rem [4j which is totally inspired by the geometry brought in our investigation by 
the theory of amoebas. Theorem [4] gives the answer to the question whether an 
average energy of an ensemble permits the thermodynamical limit. Namely, the 
domain of admissible average energies coincides with the interior of the convex 
hull of the spectrum. 

2. Amoebas of complex hypersurfaces 

For convenience we shall denote by T n the set (C \ {0})™. 

Definition 1 ([lj). The amoeba Ay of a complex algebraic hypersurface 

V = {z G T" : Q(z) = 0} 

(or of the polynomial Q) is the image of V under the mapping Log:T" — > M. n , 
determined by the formula 

Log: (zi,...,z n ) (log|zi|,...,log|z„|). 
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The term amoeba is motivated by the specific appearance of Ay in the case 
n = 2. It has a shape with thin tentacles going off to infinity (see Fig. [lj. The 
complement K™\^4y consists of a finite number of connected components, which 
are open and convex pQ . The basic results on amoebas of algebraic hypersurfaces 
can be found in HEEIUS] " 

We denote by Mq the Newton polytope of the polynomial Q, that is, the convex 
hull in R" of all the exponents of the monomials occurring in the polynomial Q. 
For each integer point v e Mq we define the dual cone C„ to the polytope Mq 
at the point v to be the set 

C v = {se K" : (s, v) = max (s,a)}. 

oc£Nq 

We recall that the recession cone of a convex set E C K ra is the largest cone, 
which after a suitable translation is contained in E. The connection between the 
combinatorics of the Newton polytope Mq of the polynomial Q and the structure 
of the complement of the amoeba Ay is described by the following result. 

Theorem QllJ). On the set {E} of connected components of the complement 
M™ \ Ay there exists an infective order function 

v. {E} -> Z" n Mq 

such that the dual cone Cw^n to the Newton polytope at the point v(E) is equal 
to the recession cone of the component E. 

This means that the connected components of the complement R n \ Av can 
be labelled as E v by means of the integer vectors v = v(E) e Mq (see Fig. [TJ. 

The value v(E) of the order function allows for two interpretations. On the 
one hand, v(E) is the gradient of the restriction to E of the Ronkin function 
for the polynomial Q (see |12|). The Ronkin function is a multidimensional 
analogue of Jensen's function and finds numerous applications in the theory of 
value distribution of meromorphic functions. On the other hand, components of 
the vector v{E) are the linking numbers of the basis loops in the torus Log -1 (a;), 
for any x € E, and the hypersurface V (see [TT] or [5]). 



1 



Mikael Passarc 



, Dmitry Pochekutov, August Tsikh 



Remark. The set vert A/q of vertices of the polytope A/q belongs to the image 
of the order function v . In other words, for each vertex f3 € A/q there is a 
component Ep with recession cone Cp (U i\14V - The existence of components E v 
corresponding to other integer points v £ A/q \ vert A/q depends on the coefficients 
of the polynomial Q. 

There is a bijective correspondence between the connected components {E„} 
of the complement W 1 \ Ay and the Laurent expansions (centered at the origin) 
of an irreducible rational fraction F(z) = P(z)/Q(z) (see [3 Sect. 6.1]). The 
sets Log~ 1 (E lJ ) are the domains of convergence for the corresponding Laurent 
expansions. One may therefore label such an expansion using the components of 
the amoeba complement, or using the integer points in the Newton polytope. For 
instance, the Taylor expansion of a function that is holomorphic at the origin will 
always correspond to the vertex of the Newton polytope A/q with coordinates 
(0,...,0). 

In Sects. [5][7] we shall see that, when working with partition functions, one 
needs to consider amoebas also of non-algebraic complex hypersurfaces. Let Q 
be a Laurent series in the variables z = (z±, ... , z n ): 

Q(z) = J2 a ^ a - 

a£AcZ n 

We assume that its domain of convergence is non-empty, and that Q(z) ^ 0. We 
shall also make the assumption that Q does have zeros in GflF. Let 

V = {zeGDT n :Q(z)=0} 

be the hypersurface given by the zeros of the analytic function Q(z). The amoeba 
for V is defined as in the algebraic case: Ay — Log(V). 

We introduce the notation Q = Log(G) for the image of the convergence 
domain G of the series Q. It is well known that Q is a convex domain. In the 
algebraic case, when Q is a polynomial, the set Q is all of W 1 , and the amoeba 
Av is a proper subset of Q . In the general case it may well happen that there 
is an equality Ay — Q. To avoid this situation, we require that the summation 
support A of the series Q lies in some acute cone, that is, the closure Af of the 
convex hull c\i(A) does not contain any lines. 

Theorem 1. If for the series Q the set Af = ch(A) does not contain any lines, 
then the complement Q \ Ay is non-empty. To the set {v} of vertices of the 
polyhedron AA there corresponds a family {E u } of pairwise distinct connected 
components of the complement Q \ Ay . The dual cone C v to M at the vertex v 
coincides with the recession cone for E v . 

Proof. Assumption of the theorem implies that the set of the vertices vert (Af) 
is nonempty. The argument is similar to the one for the algebraic case (when Q 
is a polynomial and Q = R n ) that is given in [TT] and [T3]. First one shows that 
for each vertex v £ Af a suitable translate of the cone C„ is disjoint from Ay, so 
that one can associate with the vertex v the component E v of the complement 
Q \ Ay, which contains this translated cone. Here the only difference is that, 
when Q 7^ W 1 , one must show that the translated cones are contained in Q. This 
follows from the fact that the dual cones C v at the vertices of Af all lie in the 
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cone — C V (A/"), where C v (Af) is the dual cone of the recession cone C(Af) of Af, 
together with the multidimensional Abel lemma |15) . which says that the cone 
— C V (A/") lies in the recession cone of the domain Q. 

Next, just as in .1,4], one associates to the collection of n-cycles F v = Log~ 1 (x„), 
with the point x v taken in the translate of C v , a collection of de Rham dual n- 
forms w M which are meromorphic in GflT" with poles on V . Namely, we choose 

' - ^ - dZl A...A^, M evert(A0 



(2m) n Q{z) z x z n 

(recall, that is the Laurent coefficient of Q). For points z € E v we have 
\a v z u \ > |g„(z)|, where g v {z) = Q{z) — a v z v . Hence, the meromorphic function 
l/Q(z) can be developed into a geometric progression 

1 _£ ( _ 1)fc sSW 



uniformly converging on r v , and one has 

(-i) fc r _ 

(27ri)™ Jp u a v z v \ a v z v j z\ z. 



„ M _ y' (- 1 ) fc /" f g*(*A fc ^1 A A ffcn 



r " fc=0 



The leading term of Q(z) with respect to the orders, defined by weight vectors 
from C„, is equal to a v z v . This yields that all the integrals in the sum vanish 
for v 7^ /Lt; and if v = fi, the only one nonzero summand occurs for k = and 
equals 1. Therefore, 



w p = Sup, 



and by the de Rham duality [16] the cycles r u , v g vert A/ - are linearly indepen- 
dent in the homology group H n ((G n T n ) \ V"). The cycles Log _1 (x) for x from 
the same connected component of Q \ Ay are homologically equivalent, this im- 
plies that the connected components {-Ei/}^ eV ert(A r ) are pairwise distinct. Since 
the n-dimensional cones of a fan dual to C{N) coincide with the cones C„ and 
C v C E v , one has that C„ coincides with the recession cone for E v . □ 



3. The amoeba contour and the logarithmic Gauss mapping 

In Sect. [2] we saw that certain information about the position of the amoeba of a 
complex hypersurface is given by the combinatorics of the integer points of the 
Newton polytope (or polyhedron) of the polynomial (or series) that defines this 
hypersurface. Here we shall describe an object associated with the amoeba, that 
reflects the differential geometry of the hypersurface. The study of this object 
can be carried out with more analytic methods. 

The contour Cy of the amoeba Av is defined (see [13]) as the set of critical 
values of the mapping Log: V — > R" , that is, the mapping Log restricted to the 
hypersurface V. We observe that the boundary dAy is included in the contour 
Cy 1 but the inverse inclusion does not hold in general. Note, a contour of an 
amoeba for Harnack's curve coincides with a boundary of the amoeba 0[5] (there 
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Fig. 2. The amoebas and their contours for the graphs of polynomials 1 — 2z — 3z 2 (left) and 
1 + z + z 2 + z 3 (right, the normal line I to <9i?o,i with a directional vector q and points x, y 
illustrate the proof of the Theorem 2). 



is the amoeba of the Harnack curve on Fig. |Tj) . Herewith, a real section FnK 2 of 
Harnack's curve consists of fold critical points of the projection Log : V >->• Ay. 
Fig. [2] depicts the amoebas of the complex curves, which contours do not coincide 
with the boundaries, besides that the points a, b, c and d are the images of 
Whitney pleats. 

We recall (see [2lll7j) that the logarithmic Gauss mapping of a complex hy- 
persurface V C T n is defined to be the mapping 

7 = 7vr:regV->CP n _i, 

which to each regular point z € reg V associates the complex normal direction 
to the image \og(V) at the point log(z). (Here log, in contrast to Log, denotes 
the full complex coordinatcwise logarithm.) The image 7(2) does not depend on 
the choice of branch of log and it is given in coordinates by the explicit formula 
0: 

7(2) = (ziQ' Zl (z) : ... : z n Q' Zn {z)) . 

The connection between the contour Cy and the logarithmic Gauss mapping 
is given as follows. 

Proposition 1 (|2j). The contour Cy is expressed by the identity 

Cy = Log (7^ (»>„_!)) . 

In other words, the mapping 7 sends the critical points z of Log|y to real direc- 
tion 7(2) which is orthogonal to the contour Cy at Logz. 

The inverse z — r )~ l (q) of the logarithmic Gauss mapping is given by the solu- 
tions to the system of equations 

fQ(*) = o, m 

\q n z j Q' Zj -q j z n Q' Zn =Q, j = 1, . . . , n - 1 . { ' 

For a fixed vector q € Z™ = Z n \ {0} the solutions to the system (Til) consist 
of the points z g V at which the Jacobian of the mapping (Q(z), z q jhas rank 
< 1, which means that the following statement holds. 
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Proposition 2. A point w € reg V is a critical point for the monomial function 
z q \ v if and only if the logarithmic Gauss mapping takes the value q at w, that 
is, 7(w) = q. 

Notice that if V is the graph of a function of n variables z = (zi, . . . , z n ), so 
that it is the zero set of the function Q(z,w) = w — f(z), then the logarithmic 
Gauss mapping is given in the affine coordinates Sj = qj/q n +i, j = 1, .. . ,n of 
CP„ by the formula 

f'z 

Zj^r = -Sj, j = l,...,n. (2) 



4. Asymptotics of Laurent coefficients 

Let E be a connected component of the amoeba complement with smooth bound- 
ary dE. The cone generated by the outward normals to dE will be called the 
component cone of E and denoted by Kg- It is clear that Kg is a cone over the 
image of dE under the ordinary Gauss mapping a : dE — > S n ~ l . 

Definition 2. The smooth boundary dE of a connected component E is said to 
be simple if for each x € dE the real torus Log -1 (a;) intersects V in a unique 
point, and if moreover the logarithmic Gauss mapping 7 of the hypersurface V 
is locally invertible at this intersection point. 

The following result is a consequence of Lemmas 1 and 2 in the paper (TS] , 
which exhibits a class of simple boundaries in the case where V = i~/ is the 
graph over the convergence domain of a power series f{z) = ^2 uj a z a . 

a£AcN" 

Proposition 3. If € A, the coefficients ui a are positive, and the set A gener- 
ates the lattice Z™ as a group, then the boundary of the component Eq 1 of the 
complement of the amoeba Ar f is simple. 

As it follows from an example of a polynomial / = \ — 2z\~Zz\ (see Fig. [2| the 
condition of coefficients uj a to be positive is essential in Proposition [3j Namely, 
the preimage of the inner point of the arc (a, b) C 9-Eo,i consists of two points 
on the graph /y, and the boundary point a or & have one preimage on Jy, but 
the logarithmic Gauss mapping has no inverse at a and b. 

In view of the convexity and smoothness of dE each point x € dE is the 
preimage x = a^ 1 (q) of a point q g Kb- 

We consider the expansion of the meromorphic function F = P(z)jQ(z) in a 
Laurent series 

F(z) = J2 c - za > ( 3 ) 

that converges in the preimage Log -1 (-E) of a complement component E of the 
amoeba of the polar hypersurface V — {z : Q(z) = 0} of F. For a fixed q £ 1™ 
we define the diagonal sequence c q .k = C( gi q y k of Laurent coefficients c a from 
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Theorem 2. Let the boundary dE be simple. Then for each q € Z™ n Ke the 
diagonal sequence {c q .k} has the asymptotics 

c q . k = k^ ■ z-"- k (q) ■ {C(q) + Oik- 1 )} , (4) 

as k — > +00. Here z(q) — V (1 Log- 1 (a~ 1 (q)), and the constant C(q) vanishes 
only when P(z(q)) = 0. 

Proof. The idea of the proof is to choose the cycle of integration Log -1 (a;) in 
the Cauchy formula 

1 f F(z) dzi dz n 

(2tTz)™ JLog" 1 ^) zq Zl Zn 

for those x that lie near the point y = Logz(q) g dE on the line I = {y + qt : t E 
R}, which is tranversal to dE (see Fig. [2]). In view of the assumed simplicity of 
dE, the torus Log~ 1 (y) c Log -1 (7) intersects V in a unique point, and Log^ 1 (/) 
intersects V in a neighborhood of z(q) along a (n— l)-dimcnsional chain h C V. 
By means of residue theory one shows (see [19] for the case n = 2) that, as 
a function of the parameter fc, the integral is asymptotically equivalent, as 
k — > +00, to the oscillatory integral 

2m f reswe-*W° ez > 



where 



P(z) dzi dz 7 , 
A ... A ' 



{2m) n Q{z) zi "' z„ 



and rescj = Qu/dQ denotes the residue form for uj. The phase ip — (q, log z) = 
log 2 9 has the unique critical point z(q) on h, at which Hetp attains its mini- 
mal value. A direct computation shows that the Hessian Hess tp vanishes on V 
simultaneously with the Jacobian of the logarithmic Gauss mapping. Since dE 
is simple, this Jacobian is not equal to zero at z{q), and hence z(q) is a Morse 
critical point for the phase ip. Using the principle of stationary phase we obtain 
formula Q with the constant C(q) being the value at the point z(q) of the 
function Pjz\ ■ . . . ■ z n ■ Q' Zn ■ (Hess (p) 1 / 2 . □ 



5. The thermodynamic ensemble and its most probable distribution 

We consider a thermodynamic ensemble il, consisting of N copies of some phys- 
ical system. Usually (see for instance [5], [IH]) HO]) or [|) the system is 
characterized by energy values from a spectrum 

= £ < £1 < £2 < ■ ■ ■ , Ej eZ . 

Each choice of energies in the systems of the ensemble defines a state of the 
ensemble. A basic question in the study of the behavior of an ensemble concerns 
the preferred states of the ensemble, as N — > 00. 
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We will consider a more general situation, where the system is characterized 
by a multidimensional quantity e k = (el, ...,£?) from a given spectrum 

© = {e*}*=B^ C N" , 

in which we for convenience shall assume that Eq — 0. Futhermore, we shall 
consider spectra from the lattice Z n , which lies in some acute cone in R™ D Z™. 
We introduce the quantity 

JV! 

W(a)=W(a ,ai,...)= , ,', , (6) 

a !ai!a 2 ! . . . 

expressing the number of different states of the ensemble, for which exactly 
of the systems is in the state with parameter value ek- We also say that is the 
£fc energy occupation number in the ensemble. It is clear that in ^ one should 
have 

X>* = tf, (7) 

fc 

^a fc e fe =£, (8) 
fc 

where £ = (£\, ...,£„) is the energy of the ensemble and the summation is over 
the index k that enumerates the elements Ek of the spectrum. The collection of 
numbers a — (a&) is said to be admissible if it satisfies conditions ([7| and QSJ) . 

By definition, the mosf probable distributions of energies among the systems 
of the ensemble (for N 1) correspond to those a that occur most frequently, 
that is, those that realize the maximum 

max W(a) 

a 

among all admisible collections a. 

When considering the problem of describing the most probable energy distri- 
butions one makes the assumption that the vector £/N = u is kept constant, 
that is, the average energy u — (iti, . . . ,u n ) of the ensemble systems is fixed. 
Under this condition, vector relation ([8| written out coordinate-wise gives n 
relations among the independent variables a&. Just as in the case of a scalar 
spectrum (n = 1, see for instance [7]), following the approach of Boltzmann, one 
uses the Langrange multiplier method to find the distributions that maximize 
W(a), which we write now W u (a) (see [TH] for details). The Lagrange multipli- 
ers /x j , that correspond to the coordinate- wise connections of vector relation Q , 
provide an important language for the solution of the assigned problem. More 
precisely, by introducing the partition function as the series 



z(fji) = z(m, ...,(j, n ) = ^2 



fc 

we obtain the fundamental thermodynamic relations: 



-V^logZ^u, a k =N- 
where is the gradient with respect to the variables fj,. 
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In order to apply methods from analytic function theory and the method 
of stationary phase, it is more convenient for us to consider other (complex) 
coordinates Zj — e~ Mi , j = 1, . . . , n. In these coordinates the partition function 
has the form 

k a£6 

Analogously, the fundamental thermodynamic relations assume the form 

Zjz) 

Let us give an interpretation of these relations by the following 
Statement 1. For N 3> 1 the occupation values ak = ak(u), computed from 



*k=N^—. (11) 



the formula (11) in the solutions z = z(u) of the system of equations (10), are 
the coordinates of the critical points for the function W u (a); in particular, the 
most probable distributions a — (ak) may be computed by means of the indicated 
formula for suitable solutions z(u). 

The comparison of formulas (10) and ([2| shows that the solutions z(u) to 
system (10) is nothing more, thantne inverse image 7 _1 (— u) of the logarithmic 
Gauss mapping 7 : Tz — > CP„ of the graph J/ of the partition function Z(z). 
However, the list of links between the mathematical notions introduced in the 
first and the second sections and the fundamental thermodynamic relations goes 
beyond this shallow observation. Another great of importance link can be figured 
out by the computation of critical values of the function W u (a). 

Since the logarithm is a smooth function, the critical points of W(a) and 
log W(a) coincide. The latter function can be written for large N with a help of 
Stirling's asymptotic formula in the form 

log W(a) = N (log N - 1) - a k (log Ofc - 1) . 

k 

The critical values of this function (under restriction £ /N — u) are 
logW u = log [z(u)- £ Z(z(u)) N ] =N (log Z(z(u))- («, log «(«))). (12) 

It is easy to check this equality by substitution in the previous expression for 
log W(a) of values ( |TT] ) for a^ evaluated at the solutions z — z(u) to the sys- 
tem (10), using relations (|7| and |8]). 



We are interested in the critical values logW„ only for real u, i.e. u £ M™. 
The portion of a critical value attributed to one system of an ensemble, i.e. the 
value ^ 

S u -■ —logW u = log Z(z(u)) - (u,logz(u)) 

plays a role of entropy. Since in the logarithmic scale log z = — fx one has 

u = - V p log Z = Vio g z log Z, 

the entropy S u considered as a function in variables u is the Legendre transform 
of the logarithm of a partition function in the logarithmic scale. 
Thus, based on Proposition [l] we get the following 
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Fig. 3. Amoebas for the graph of the partition function 1 + z 2 /(l — z): the full amoeba (on 
the left), its pieces over \z\ < 1 and over \z\ > 1 (in the middle) and with a common tangent 
segment [b, b'] to 9i?o,i and dE\ } \ (on the right). 



Statement 2. The liftings of the solutions z(u) of the system (10), for u G 
W 1 C KP„, to the graph Tz of the partition function coincide with inverse image 
7 _1 (— u) of the Gauss logarithmic map 7 : Tz ~ > CP„. On the amoeba Ar z of 
the graph these solutions parametrize the contour of the amoeba. The values S u 
of the entropy coincide with the critical values of the linear function 

lu(%*) -^n+l U\X\ ... U n X n7 

restricted to the boundary OEq A of the connected component Eqi of the comple- 
ment R n+1 \Ar z - 

For certain spectra 6 the partition function Z(z) admits an analytic contin- 
uation outside the domain of convergence of its series representation (JoJ) with 
new "twin spectra" & C Z™ appearing. Let us consider two examples. 

Example 1. The partition function Z for the spectrum 6 — {0, 2, 3, 4, . . .}, n = 1, 
is equal to the rational function l + z 2 /(l — z), which outside the unit disk \z\ < 1 
has the development 

z = -(*+i + l + . ..) = -£*■ 

aee' 



z 



where & = {1,-1,-2,...}. We can consider the thermodynamic relations 
(flO|,([TTj) also in the complement {\z\ > 1} of the unit disk. The correspond- 
ing pieces of the amoeba of the graph of this rational function are depicted on 
Fig. [3] in the middle. 

On Fig. [3] the points a and d depict the points at infinity, where the normal 
vector [—« : 1] to the contour of the amoeba at a and c' equals [0 : 1] and [—1 : 1] 
respectively. The boundaries dE Q i and dE\ \ have a common tangent at the 
points b and b' (a simple computation shows that the normal vector [— uq : 1] 
corresponds to the value uq — 1/2). The set of the normal vectors to the arc 
(a, b) C dEo.i coincides with the set of the normal vectors to (a', b') £ dEi,x, 
the same holds for the pair of arcs (6, c) and (&', c'). The tangents at the points 
of the arc (b, c) lie higher, than parallel to them tangents at the points of (&', c'), 
the tangents at the points of (a, 6) and (a' ,b'). It follows from Statement [2] that 
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the maximal value of the entropy S u for < u < uq corresponds to a solution 
z(u) projected on the arc (a', 6'), and it is from the domain {\z\ > 1}. 

However, the combinatorial interpretation of W(a) forbids us to consider the 
domain {\z\ > 1}, because all the occupation numbers in ([TT]) for z > 1 and some 
of them for z < — 1 are negative. Moreover, the partition function is negative at 
the points that project on the boundary dE^-y. 

The next example shows that in several dimensions we can overcome such 
limitations. 



Example 2. Consider the spectrum 

6 = {(0,0)}U{(2,2)+§} + {(4,4)+§}, 

where S is the semigroup (2, 1) • N + (1, 2) ■ N (see Fig. gon the left). 

The partition function X^aee z ° converges in the domain D = {\zfz2\ < 1, 
|ziz|| < 1} and equals 

Z (z) = i + d + ^fe 2 



{i-ziz 2 ){i-z lZ ir 

The development of Z(z) in the domain D' = {\z\z2\ > 1, \ziz$\ > 1} is again a 
partition function, i.e. it is a power series 

Z(z) =J2 z ° 

with summation over the spectrum 

6' = {(0, 0)} U {(-1, -1) - §} U {(1, 1) - §} 

(see Fig. |i]in the middle). 

The full amoeba of the graph Fz corresponds to the polynomial 

(w - 1)(1 - zfz 2 )(l - z x z\) - z\z\ - z\z\ 

in three variables Z\,Z2 : w. Points (0, 0, 1) and (3, 3, 1) are vertices of the Newton 
polytopc for this polynomial, therefore the complement to the full amoeba of the 
graph Fz contains connected components -Eo.o.i and ^3,3,1- Since the Laurent 
coefficients of the developments of Z in the domains D and D' are positive, the 
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boundaries 8Eq q,i and dE 3t3 ,i are the Log-images of the graph Fz over the real 
domains D n and D' n . Consider the "diagonal" function 

The amoeba of its graph can be embedded in the amoeba Ar z by the mapping 

i : (log |*|, log \Z(t, t)\) ^ (log(i), log(t), log \Z(t, t)\). 

The boundaries of the components £0,1 and Eq_i of the complement to the 
amoeba of the graph of Z(t, t) are the Log-images of pieces of the graph over the 
intervals < t < 1 and 1 < t < 00, respectively. The amoeba Ar z lives in the 
space M. 3 of variables Xi,X2,Xs; and the plane X\ = x-i cuts out in the surfaces 
9£o,o,i and 3,1 two pieces, the images i(9i?o,i) and i(dE 6 ^), respectively. 
As in ExampleJTJ the curves dEo^i and dE§^ have a common tangent line, lying 
below these curves, since they are convex. 

In view of the symmetry of Ar z with respect to the plane x\ = x%, there exists 
a common tangent plane r to surfaces dEo,o,i and dE^^ t \ with the property 
that t crosses the common tangent line to the embeddings i(dEo t i) and i(dEgi) 
symmetrically with respect to the plane xi — x 2 . As it follows from results of 
Sect.[7j the vector [u\ : 112 : 1] is normal to the tangent plane r, if u = (1*1,1*2) 
belongs to the intersection of interiors of convex hulls of the spectra & and &' , 
i.e. to the double-shaded rhombus on the right of Fig. El In general, the rhombus 
is divided by some curve 7 into two domains, such that the value of the entropy 
S u (corresponding to the ensemble with the spectrum 6) is greater than that 
of the entropy S' u (corresponding to the ensemble with the spectrum ©') in the 
first domain and is less in the second one. Perhaps, this phenomenon may be 
considered as a tunnelling transition from one ensemble to another in a way to 
increase the entropy, when we choose the value of the energy u on 7. 

At the end of this section, we show that the notion of multidimensional spec- 
trum, our starting point, leads to the notions of the multidimensional temper- 
ature and the vector of thermodynamic forms. For this purpose, we compute 
the differential of logarithm of a partition function assuming that the variables 
Zx, . . . , z n are positive and entries of the spectrum {sf.} vary in some neigh- 
bourhood of lattice points in E™, i.e. we consider the spectrum {e^} to be vari- 
able. 



In accordance with (10) and (11) 



Z' d7 Z 'l 



k j 



= (u,dlog z) + ^2Y1 ~Z l0g Z 3 dS k = ( U > Z ) + ( l0 g Z > Jt akd£k / 

j k \ k I 

Hence, we get the following expression for the differential of the entropy 
dS = d(logZ- (it, log*)) = (-log z,du) + (log z, ^^2 a kde k ^j = {^f^ UJ / 
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where 

w = (wi,...,w n ), T = (Ti, . . . ,T„) 

denote the vector of the thermodynamic forms and the vector of the temperature 
with components 



= duj - 1/N^a k de k , T = -1/log; 



6. The average value a of the admissible collections {a} 

In the preceding section we gave a description, following Botzmann, of the most 
probable distributions of the ensemble. However, the method that was used is 
somewhat limited, since the extremal points ( 11 ) for Q are obtained by applying 



the Stirling formula to a^!, and this is only justified for large values of a k . In the 
case of a scalar spectrum, the Darwin-Fowler method offers a possibility to avoid 
this drawback. It consists in a description of the asymptotics of the averages of 
the occupation numbers. We shall analogously describe the asymptotics of the 
averages of the occupation numbers, when the energy spectrum is composed of 
vector quantities. In this section we show that this problem is equivalent to the 
problem of describing the asymptotics of the diagonal coefficients of a Laurent 
expansion of the meromorphic function w/(w — Z{z)). 

Definition 3 (|9j, [7j)- The average value of the admissible collections {a} is 
the collection a = (Sfc) of numbers 

- _ Ea^W(a) 

where the summation is over all admissible collections a = (a^). 
For the study of the averages Uk we introduce the sum 

£ W ^ W) = £ a^Za^.. ^< -< k - (13) 

a a 

over all admissible collections a — (a^). Here the uij are real parameters, all 
varying in a small neighborhood of 1. We remark that W(a,I) — W(a), where 
I = (1, 1, . . .) is the all ones vector. Hence, for u) — I the quantity ( JT3] ) expresses 
the total number of states of the ensemble. It is not difficult to see that 



d 



"i, = - — log^2w(a,ui) 

a 

As in [S] and [7] one proves the integral representation 



(14) 

1 



dz 

^, (15) 
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where T r = {\z\\ = ri, . . . , \z n \ = r„}, and the rj are chosen so small that on T r 
one has convergence of the series 



i( z ) = /<x w ) = E ujkz " k = E ujkz ^ k ■ ■■■■ z n ■ 

k k 

Since f(z,I) = Z{z) we refer this series to be a variation of partition function. 
Since the condition < ujk < 1 + S is fulfilled, the domain of convergence G' of 
this series is non-empty and contains the origin z = 0. 
We now introduce the function of n + 1 variables 

F(z,w) = - "' 



«-/(*)' 

which is meromorphic in the domain G = G' x C w . The polar hypersurface of F 
is the graph 

r f ={(z,w)eG:w = f(z)}. 

Due to the fact that £o = 0\ the closure AT of the convex hull of the summa- 
tion support of the series w ~ f(z) contains the vertex v = (0,1). According 
to Theorem fll this vertex corresponds to a connected component E^ x of the 
complement of the amoeba Ay ■ Using a geometric progression we expand F in 
a Laurent series, convergent in Log _1 (i?o i i) C {(z,w) e G : \w\ > |/(z)|}: 

N=0 N £ 

For the Laurent coefficients Cs,-n of this series one has the integral repre- 
sentation 

If w -£ n h dzj dw 

where x € Eq x . Performing the integration with respect to w in this last integral, 
we immediately obtain (15 1. 

We thus find that the problem of describing the asymptotics of the sum (p~3| ) 
is equivalent to the same problem for the coefficients Cs,-n of the series (UW, 
for £ = u ■ N , with u being the vector of average energies. That is, it amounts 
to finding the asymptotics of the diagonal coefficients C[ Ui -\yx with direction 
vector q = (u, —1). 



7. The asymptotics of the average values 

Let the point (z*,u>*) on the graph Jy of the variation / of partition function 
be such that Log(z >t ,w*) € d-Eg.i- 

Since 8Eq i is a part of the amoeba contour, the first coordinates of the 
given point on the graph satisfy ([2| for some u 6 K.™, and the coordinate 
is uniquely determined by z*. As we let uj tend to the vector I = (1, 1, . . .), we 
get / — > Z, and the point (z*,w*) moves to the point (z,w) — (z(u),w(u)), 
whose logarithmic image lies on the boundary dE^ x of the component E$ j of 
the complement to the amoeba of the graph Fz = {w = Z(z )\ o f the partition 
function of the ensemble. Besides that, z(u) satisfies system (|l0]) . 
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Theorem 3. Suppose that the spectrum 6 — {sk} generates the lattice Z n , and 
that the point z — z(u) £ R™ satisfies the system {10). Then, as N — > oo, the 
average values ~a~k for the occupation numbers of energy Sk has the form 



a k ~ N 



Z(z) 



(17) 



and they coincide with most probable values of ak- 



Proof. By assumption the spectrum & generates the lattice IT 1 and hence, ac- 
cording to Proposition |3] the boundary dE^ i is simple. Therefore we can apply 
Theorem [2] for the asymptotics of the diagonal sequence of Laurent coefficients 
of the series ( fl6| : 

C ( „,-i).jv ~ C{q) ■ N-9 ■ (z- v (u)w*(u)) N , N -» +oo . 

Hence, taking into account the summary in Sect.[6j we find that the asymp- 
totics of the total number of states, as N — > +oo, has the form 



£ W(a t u) ~ C(q) ■ AT-f • (z- u (u) ■ f(z*(u))) N . 

a 

Now, direct calculation leads us to the asymptotic equality 



where ip = log(z _ll /(z)) denotes the phase (see the proof of Theorem J2J). In the 
right hand side of the last formula the first term is equal to zero, because is a 
critical point for the phase (p. Therefore, setting oj — I, we get from the formula 



(14) the desired asymptotics (17 1. □ 



Let us now raise the question about what the admissible values are for the 
vector u of average energies, that guarantee the existence of a solution z(u) G W] 



to the system (10), and hence provide the asymptotics (17). 

In the work of Darwin and Fowler [5] , [10] this question was not considered. 
Apparently, it was first paid attention to in Sect. 4.5.1], where it was observed 
that if the partition function is a polynomial of degree d, then the admissible 
average energies must be taken within the interval < u < d, that is, in the 
interior of the convex hull of the numbers = £o < e i < • • • < £ k = d. 

The raised question is answered by the following theorem, where we use the 
notation Af° for the interior of the convex hull in W 1 of the spectrum & — {£&}. 



Theorem 4. Suppose that the spectrum & = {sk} generates t he la ttice Z™. Then 
for every value of the average energy u € Af° the system has a unique 

solution z — z(u) in R™ , and hence for u G Af° the average values Ofc coincide 
with the most probable ones. 
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Proof. Lifting the solutions z(u) of the system of equations (10 1 for u G R" to 
the graph Ez of the partition function of the ensemble, we obtain the criticial 
values for the mapping Log|y. On the amoeba Ar z of the graph, these solu- 
tions parametrize its contour. In particular, the solutions z(u) € M™ that are 
of interest to us parametrize the boundary of the complement component Eq±. 
Thanks to the fact that the spectrum generates the lattice Z n , we know from 
Proposition |3| that to each point on 8Eq 1 there corresponds a unique vector 
q € Ke s 1 - Therefore, in order to obtain all solutions z(u) from K™ one must go 
through all vectors q from the component cone Eq ±. 

By Theorem [T] the recession cone of the component Eq x is the dual cone to 

Af at the vertex v = (0, 1), where N denotes the closure of the convex hull of 
the summation support of the series Q — w — Z(z). (See Figure ^1 where the 
recession cone is bounded by dashed lines.) The outward normals ofthose facets 
of the polyhedron N that come together at the vertex v span this dual cone. 
Therefore, the sought cone Ke 5 x is spanned by the edges of N that emanate 
from the vertex v, and thus Ke 6 1 consists of all vectors of the form q = (u, —1), 
with u G W°. □ 



We conclude with some remarks and illustrations to Theorem |4j First, the 
statement of the theorem still holds if one shifts the spectrum 6 by a noninteger 
vector. For example, the domain of admissible average values of energy in the 
case of the Plank oscillator with the spectrum {1/2 + N} equals {u > 1/2}. Such 
domain for the Fermi oscillator with the spectrum 6 = {0, 1} is the interval 
{0 < u < 1} (see pj ch. 4]). The latter case is depicted on the right of Fig. [5] 
Example[2]of Sect.[5]deals with the "twin-spectra", and the sectors on Fig.[3]are 
the domains of admissible average values of energy in the corresponding cases. 
These sectors have a nonempty intersection, the double-shaded rhombus (Fig.[2J 
on the right). 
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